Acoustic Plasmons in MgB 2 



O 

o 



K. Voelker 1 , V. I. Anisimov 1 ' 2 ) and T. M. Rice 1 
1 Theoretische Physik, ETH Honggerberg, CH-8093 Zurich, Switzerland 
2 Institute of Metal Physics, Russian Academy of Sciences, Ekaterinburg, Russia 

(February 1, 2008) 

We present strong evidence for the existence of an acoustic plasmon mode, that is, a quadrupo- 
lar charge collective mode with linear dispersion, in MgB2. This mode may be responsible for the 
anomalously small value of the Coulomb pseudopotential required to explain the high superconduct- 
ing transition temperature of 40K. 
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The recent discovery of superconductivity in MgB2 [1] 
has spurred a wave of activity both on the experimen- 
tal and on the theoretical side. The transition tem- 
perature of T c — 39. 5K is by far the highest found in 
any binary compound. With the exception of Cgo, it 
is surpassed only by the high-temperature copper-oxide 
superconductors. A detailed analysis of the electron- 
phonon interaction finds relatively strong coupling to a 
high frequency phonon mode [2]. However, the observed 
T c = 39. 5K could be obtained only by assuming that 
the Coulomb pseudopotential takes an anomalously small 
value, \i* ~ 0.02 — a value roughly seven times smaller 
than for usual s-p superconductors. Here, we propose 
that an acoustic plasmon mode exists in MgB2. This 
leads to an attractive contribution, thereby lowering the 
value of the Coulomb pseudopotential. 

The idea of acoustic plasmons goes back to Pines [3], 
and, following Frdhlich, has been considered as a mech- 
anism for superconductivity in the transition metals in 
the past [4]. In the presence of two carrier species with 
very different effective masses, the light carriers can act 
to screen the Coulomb repulsion between the heavy car- 
riers, and thereby lower the plasma frequency of the 
heavy carriers. The resulting plasmon mode will have 
a linear dispersion uj a — v a q (hence the term acoustic 
plasmon), where v a is of the order of (but somewhat 
higher than) the Fermi velocity of the heavy carriers. 
The acoustic plasmon mode is Landau damped by decay 
into the electron-hole continuum of the light carriers, but 
for small wave vectors q the mode will be located within 
the lower range of the electron-hole continuum, so that 
the damping may be weak, and the acoustic plasmon may 
exist as a well-defined collective mode. 

Since no experimental evidence of acoustic plasmons 
could be found, and one was able to explain superconduc- 
tivity in the transition metals within the usual phonon 
mechanism, interest in the topic began to wane in the 
late 1980's. We believe that MgB 2 is a good candidate 
for further research on acoustic plasmons and their con- 
tribution to superconductivity, due to its Fermi surface 
characteristics. These have been first reported in Ref. [5], 
and will be briefly described below: 

MgB2 consists of honeycomb layers of boron atoms 




FIG. 1. Fermi surface of MgB2. The figure is taken from 
Ref. [5]. Holes in the cr-band form cylinders around the 
rA-line. The 7r-band has electron and hole pockets located 
near the H- and -ff-points, respectively. 



that are stacked vertically with no displacement. The 
magnesium atoms are located at the center of the 
hexagons formed by boron, but between the boron 
planes. The resulting structure is called the AIB2 struc- 
ture. The in-plane lattice constant is a — 3.086A, while 
the interlayer spacing is c = 3.524A. Hence there is a 
strong anisotropy in the boron-boron bond lengths. To 
a good approximation, the outer shells of the Mg atoms 
can be assumed as being fully ionized, so that the boron 
system bears a certain resemblance with graphite, both 
structurally and electronically. The differences in the 
electronic properties of MgB 2 and graphite are mainly 
due to the attractive potential of the Mg ions, which en- 
hances the dispersion of the p 2 -orbitals parallel to the 
c-axis through hybridization with the Mg s-orbital. 

Carriers are situated in two hole bands derived from 
the tr-bonding boron p^^-orbitals, that are essentially 
two-dimensional, and in one electron and one hole band 
derived from the 7r-bonding boron p^-orbitals (see Fig. 1). 
The smaller overlap of the latter in the x, y-directions 
compensates for the larger interatomic spacing in the z- 
direction, so that the 7r-bands are nearly isotropic. 

An essential feature of the band structure is the large 
difference in the z-components of the Fermi velocities be- 
tween the a hole bands and the tt electron and hole bands. 
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As will be shown explicitly below, the existence of two dif- 
ferent velocity scales gives rise to a weakly damped acous- 
tic plasmon mode that corresponds to long-wavelength 
charge density fluctuations between these bands. While 
the total charge density p n + p a does not fluctuate, the 
charge difference p T — p„ will vary locally. Since the 
cr-orbitals are mainly located within the Boron planes, 
while the charge in the 7r-orbitals is situated between the 
planes, this charge transfer between the orbitals gives rise 
to quadrupolar charge fluctuations. 

We will now show the existence of a well-defined acous- 
tic plasmon mode within the random phase approxima- 
tion (RPA) . For our purpose it is convenient to define the 
polarization insertion Tl ab for different carrier species a 
and b as indicated in Fig. 2. For a single carrier species, 
n aa then reduces to the standard definition [7]. The po- 
larization insertion is related to the time-ordered density- 
density correlation function C ab between different types 
of carriers as 

C ab (r,t) =-i(Tn,(r l i)ni(0 1 0)> = hU ab (r,t). (1) 

Here q — (q, uS) is a momentum four- vector, Uo(q) = 
Ane /q 2 is the bare Coulomb potential, and n a — n a — 
(n a ) denotes the fluctuations in the charge density of 
carrier species a. Specializing to the case of two carrier 
species of opposite charge (electrons and holes) , we can 
express Tl ab in the random phase approximation as 



8 ab IL° a (q) - U (q)Tl<l( q )n°( q ) 



(2) 



1 - U (q)Ul(q) - U (q)H° 2 (q) 

where the lowest-order polarization insertion 11^ is given 
in terms of the non-interacting Green's functions as 

2 f d^k 



ih J (2tt) 



jG^MGHk + q). 



(3) 



We now define the total charge fluctuations n_ = fi\ — 
n-2 and the charge transfer fluctuations n + = fi\ + n-z. 
Then the corresponding correlation functions are given 
by 



lr m = n°(g) + n°(g) 

h W l-U (q)Il<i(q)-U (q)Il»(q)> 



(4) 



and 



i r M _ ng(g) + n°( g ) - 4t/ (g)n;(g)n°(g) 
h +[q> i-EM9)n2( 9 )-tfo(«)nS(?) ' 



njg) = 




FIG. 2. Definition of the polarization insertion for different 
carrier species a and 6 



0.03 r 



0.02- 



0.01 - 



q = (0, 0, 0.01 Jt/c) 



s 



-0.01 



-0.02 



-0.03 




FIG. 3. Lowest-order polarization insertion II as a func- 
tion of frequency, for a fixed wave vector in the z-direction. 
The black solid and dashed lines show the real and imagi- 
nary part of the total lowest-order polarization. The black 
arrow indicates the frequency of the acoustic plasmon mode, 
at which ReII (aj) vanishes. Also shown are the individual 
contributions from the two a hole bands (<7i, blue and 02, 
green), and from the 7r electron (ni, red) and hole (II2, pur- 
ple) bands. Solid lines indicate the real parts, and dashed 
lines the imaginary parts. Units for q and oj are A and eV, 
respectively. 



A collective mode corresponds to a pole in the corre- 
lation function in the complex frequency plane, and can 
thus be identified by the condition 
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(6) 



which remains valid for more than two carrier species 
as well. For small wave vectors Uo(g) will be very 
large, so that we can set the left-hand side to zero in- 
stead. Inserting the above condition into the numerator 
of Eq. (4), we see that the pole's residue is proportional 
to l/U~o(q) ~ q 2 , so that fluctuations in the total charge 
density are suppressed for small wavevectors. Due to the 
additional term in the numerator of Eq. (5) this is not 
the case for the charge transfer fluctuations described by 
C+. 

We calculated the band structure with the tight- 
binding linear muffin-tin orbital (TBLMTO) method [8], 
and used the following analytical fits for our calculations 
of I! : For the a bands, 



e(k) 



24 ct) cos(fc z c) 



2m 



1.2 



-{kl + kl), (7) 



where e 







0.66 eV, t 



0.22 eV, mi = 0.27 m e , and 



ni2 = 0.50 m e , with m e being the bare electron mass. For 
the 7T bands we used 
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e(k) = 4" } - 2t< 1 ° cos(fc z c) 



(8) 



'31 1 
±2t x \j - + - cos(k x a) + cos(-k x a) cos( — fc y a), 



where = 0.16 eV, t { ? ] = 1.12 eV, t x = 1.73 eV, and 
the positive or negative sign corresponds to the electron 
and hole band, respectively. 

Within the approximate dispersion given by Eq. (7), 
the lowest-order polarization insertion of the a bands can 
be calculated exactly for the special case where q is par- 
allel to the 0-axis: 
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(qz,w) = 



mi, 2 j 6{l-Y 2 )-i6{Y 2 - I) 



nfi 2 c 



- 1 



(9) 



where Y — 4i£^ sia{qc/2)/%u, and 9{x) is the Heaviside 
step function. For q off the z-axis and for the 7r-bands 
we resort to numerical calculations. Results are shown in 
Fig. 3. From the condition (6) we can identify an acoustic 
plasmon mode with dispersion uj a = v a q, where the plas- 
mon velocity is v a ~ 1.78 AcV = 2.70 x 10 7 cm/s. The 
acoustic plasmon mode is weakly damped by decay into 
the 7r electron-hole continuum. The remaining two zeroes 
of Re I! lie deep within the electron-hole continuum and 
correspond to ovcrdamped modes, as can be seen from 
the large imaginary part of II . To estimate the damping 
coefficient A of the weakly damped mode we assume that 
the contribution from the 7r carriers to the real and imag- 
inary parts of IT is constant in the vicinity of u) a . Then 
the imaginary part of the pole's location in the complex 
plane can easily be calculated from the analytic expres- 
sion (9), and yields a Q-factor Q = o> a /A ~ 15. The 
acoustic plasmon mode is hence well-defined. 

The above analysis leads us to propose that a well de- 
fined acoustic plasmon mode exists for small values of 
the transverse component qj_. Physically this mode cor- 
responds to quadrupolar charge fluctuations involving a 
local charge transfer between a- and 7r-orbitals of the 
boron atoms. This collective mode in turn leads to an 
attractive retarded contribution to the Coulomb interac- 
tion in a window in (q, w)-space, which will modify the 
repulsive Coulomb pseudopotcntial. Further analysis is 
required to determine if the modification is strong enough 
to give the proposed reduction to a value of [i* = 0.02. 
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